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'^ , We validate a new law for the aspect ratio a = H/L of vortices in a rotating, strati- 

fied flow, where H and L are the vertical half-height and horizontal length scale of the 
^^ I vortices. The aspect ratio depends not only on the Coriolis parameter / and buoyancy 

• ■ (or Brunt- Vaisala) frequency N of the background flow, but also on the buoyancy fre- 

fj I quency Nc within the vortex and on the Rossby number Ro of the vortex such that 

a = f [Ro (1 + Ro)/{N^ — 7V^)]^/^. This law for a is obeyed precisely by the exact equi- 
^»Lj| librium solution of the inviscid Boussinesq equations that we show to be a useful model 

JS^ ■ of our laboratory vortices. The law is valid for both cyclones and anticyclones. Our anti- 

M^l cyclones are generated by injecting fluid into a rotating tank filled with linearly-stratified 

salt water. The vortices are far from the top and bottom boundaries of the tank, so there 
is no Ekman circulation. In one set of experiments, the vortices viscously decay, but as 
^ ■ they do, they continue to obey our law for a, which decreases over time. In a second set 

of experiments, the vortices are sustained by a slow continuous injection after they form, 
so they evolve more slowly and have larger |i?o|, but they also obey our law for a. The 
law for a is not only validated by our experiments, but is also shown to be consistent 
^^ with observations of the aspect ratios of Atlantic meddles and Jupiter's Great Red Spot 

^—s \ and Oval BA. The relationsh i p for a is derived and examined numerically in a companion 

^Sj ■ paper bv iHassanzadeh et al\ (J2012h . 
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1. Introduction: vortices in stratified rotating fiows 

The Great Red Spot and other large vortices such as the Oval BA (JMarcuslll993l) 



are persistent giant anticyclonic vortices in Jupiter's atmosphere. Their vertical aspect 
ratios a = H/L lie in the range 0.03 ^ a ^ 0.1. In the Atlantic Ocean, meddles are also 
long-lived anticyclones made of water of Mediterranean origin that is warmer and saltier 
than the ambient Atlantic. Their lifetimes can be as long as several years, and they have 
a ~ 0.01 (with H ~ 0.5 km and L ~ 50 km). Presumably, the aspect ratios of these 
vortices are the result of a competition between rotation and stratification. A rapidly 
rotating flow, parameterized by a large Coriolis parameter / and small characteristic 
azimuthal velocity Vg (i.e., small Rossby number Ro = Vg/fL) is controlled by the 
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Taylor-Proudman theorem: flows have httle variation along the rotation axis and form 
columnar vortices with large a. In contrast, strongly stratified flows, parameterized by 
large buoyancy frequencies 



inhibit v ertical motions and form baroclinic vortices, often appearing as thin "pancake" 
vortices ( Billant fc ChomaabOOl ) . Here p is the fluid density, g is gravity, z is the vertical 



coordinate. 

The goal of this paper is to investigate this competition and determine and verify a 
quantitative law for a. Previous experi mental studies that used constant-density fluids 
to simulate meddles or Jovian vortices (jSommeria et aZj|l988t lAntipov et al\\l98w pro- 



hibited this competition and resulted in laboratory vortices that were barotropic Taylor 
columns that extended from the bottom to the top of the tank. Thus, a was imposed by 
the boundaries of the tank. In our experiments, the vortices are created near the center 
of a large tank so that the tank's boundaries have no or little influence on a, and there 
is little or no Ekman circulation in the tank and no Ekman spin-down of our vortices. 



2. Aspect ratio law of a model vortex 

The dissipationless, Boussinesq equations for a fluid with mean density pa in the ro- 
tating frame (and ignoring the centrifugal force) are 

^ = -(v.V)v-^ + /vxz-(^^M,z and | = -(v.V)(p-p-) + .. po - (2.1) 
dt Pa Pa dt g 

where v is the divergence-free velocity, Vz the vertical component of the velocity, p is 
the pressure, z is vertical coordinate, z is a unit vector, and an overbar above a quantity 
indicates that the quantity is the equilibrium value in the undisturbed, linearly-stratified 
(i.e., with constant N) fluid. The unperturbed solution has v = 0; p = —{N'^/g)pa z + po; 
and p = (iV^po/2) z^ — g po z +pa, where pa is an arbitrary constant. One steady solution 
of the Boussinesq equations that consists of an isolated, compact vortex with solid body 
rotation il has v = v — 0; p = p; and p = p everywhere outside the vortex boundary; 
while inside the vortex: vg = fir; Vz = Wr — 0', P = —{Nc /g)poz + po; and p = 
{N^ Po/2) z'^—g Pa z+[Vl (fl + f) po/2] r'^+Pc, where r is the cylindrical radial coordinate, 
Vr and Vg are the radial and azimuthal components of the velocity; Pc is a constant equal to 
Pq~Po ^ (^+/) ^^/2, and L is the radius of the vortex at z = 0. This vortex has a uniform 
buoyancy frequency Nc throughout the entire vortex (in general, the subscript c means 
the quantity is to be evaluated at the vortex center, i.e., at r = and z = 0). The vortex 
boundary is determined by requiring that the pressure be continuous throughout the 
flow. But note that the density and velocity are discontinuous at the vortex boundary). 
Continuity of p requires the vortex boundary to be ellipsoidal: 

{r/Lr + iz/Hr^l, (2.2) 



where the semi-height is H = w2(po —Pc)/[po {Nc — N^)], and the semi-diameter is 
L = y/2 {po — pc)/['^po (/ + ^)]. This vortex has aspect ratio 
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where the Rossby number Ro = uJc/2f = fl/ f- We have defined Ro in terms of the 
vertical vorticity uj at the vortex center to make our expression in equation (|2.3|) con- 
sistent with a more gener al relationship for a that is derived in the companion paper of 
Hassanzadeh et al\ ( 20121 ) by assuming that the vortex has cyclo-geostrophic balance in 
the horizontal directions and hydrostatic balance in the vertical direction. For cyclones: 
Ro > 0, Po > Pc and Nc > N. For cyclostrophic anticyclones: Ro < — 1, po > Pc, and 
Nc > N, but for quasi-geostrophic anticyclones: — 1 < Ro < 0, pq < p^, and Nc < N. 



3. Comparison with previous observations or predictions of vortex 
aspect ratio a 

Our law (|2.3p for the aspect r atio a agrees w i th pre vious experiments in rotating 
and stratified flows. For example, iBush fc Woods! ( 19991 ) created vortices in a rotating 
stratified fluid from the break-up of small-diameter rising plumes. Very little ambient 
fluid becomes entrained in a rising plume or when the plume breaks and rolls up into 
vortices, so the vortex cores have nearly uniform density, so Nc = 0. Due to angular 
momentum conservation, the vorticity within plumes rising from a small diameter orifice 
have no angular momentum or angular velocity when viewed in the incrtial frame. When 
viewed in the rotating frame with angular velocity //2, these plumes have an angular 
velocity ft — — //2. Thus, in the rotating frame these vortices are anticyclones with 
Ro ~ —0.5. With parameter values Nc = and Ro ~ —0.5, equation (12.311 pr edicts 
a = 0.5f/N, which agrees well with the experiments of iBush fc Woodsl ()l999( ) that 



found a = 0A7f/N. 

It i s often claimed that quasi-geostrophic (QG) vortices obey t he scaling law H/L = 
f/N (JReinaud et aLl l2003':'Dri tschel et aLlll999l : lMcWimamslll985l) . independently of the 
values of Ro and Nc- This scaling is broadly used among the oceanographic community 
for vortices created from noise but is sometimes applied to Atlantic meddles and other 
vortices without close neighbors. We show in sections 4 and 5 that the aspect ratios of 
our experimental vortices and meddles have shapes that strongly depend on Ro and their 
internal stratifications Nc, and the refore a ^ f/N for these vortices. 

In a theoretical inviscid analysis, [GilJ ( 19811 ) proposed that a should be proportional to 
Ro{f /N). In this study, the velocity and the density fields were searched as continuous 
solutions of the thermal wind equation in 2D. A direct consequence of this constraint 
(that we avoid in our analysis as discontinuous fields are solutions of the primary non 
viscous equations) is that the scales of the vortices are in fact imposed by the shape of 
the solutions that match the fields inside and outside the vortices. These scales are not 
the characteristic lengths of a real discontinuous patch of vorticity and density anomaly 
as defined in section 2. Note that in reality, viscosity will smooth the discontinuity of 
the patch, but this process will no t influence the aspect ratio of the vortex. Moreover, as 
shown in the companion paper of Hass anzadeh et all ( 2012^ . Gill's law is only valid for 
the outer field contrary to his solution for the inner field that sa t isfies our law (|2.3p in 
the geostrophic limit. There fore, we thi nk that iHedstrom fc Armil (|1988l ) experimentally 
checked the law derived bv iGilll ( 1981 ) for the outer field only. It appears that in this 
experimental work, the measurements were made while the vortices were very young, 
still undergoing geostrophic adjustment and far from equilibrium. As described below, 
by performing similar experiments but waiting until the vortices come to a quasi-static 
equilibrium, the aspect ratio of the vortices indeed obeys equation (j2.3p . 

Note finally that our scaling law, modified for use with discrete layers of fluid rather 
than a continuous stratific ation, al s o app lies to the models of anticyclonic ocean eddies 
used bylNol (|l98lh and bv lGartonl (l200lh . 
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Figure 1. Aspect ratio and Rossby number. Upper panels: Determination of the aspect ratio 
Q = H/L. (a) Side view of laboratory anticyclone with iV = 2.3 rad/s and / = 2 rad/s. (b) Image 
processing of vortex in panel (a) to determine a. Here a — 0.3 with H = 1.9 cm and L — 6.3 cm. 
Lower panels: Determination of Ro. (c) Top view of laboratory anticyclone with N = 2 rad/s, 

/ = 2 rad/s, and Ro = —0.13. (d) Azimuthal ve (— ) and radial Vr ( ) velocities of flow in 

panel (c) as functions of radius r. The core has a solid body-like rotation. The magnitude of 
Vr is consistent with the R.M.S. fractional uncertainty of vg, which is ±5-10%. The fractional 
uncertainties in N are ±10%. 



4. Application to Laboratory Vortices 

We have carried out a study of vortices in a rotating, stratified laboratory flow in a 
transparent tank of dimensions 50 x 50 x 70 cm ^ mounted on a rotating table. Using 
the classic double-bucket method with salt water (|Osteiill965[ ). each experiment initially 
has a linear vertical stratification, with a constant buoyancy frequency N independent of 
location. The values of iV in our experiments varied from 1 rad/s to 2.3 rad/s. The Coriolis 
parameter / can be as large as 7 rad/ s. Two sets of experiment s ha ve been carried out 
follow ing respectively seminal works bv I Griffiths fc LindenI ( 19811) andlHedstrom fc Armil 
( 19881 ). In the first one, once the fluid in the tank reaches solid-body rotation, we briefly 
inject a small volume of fluid with constant density po through a 3 mm diameter pipe 
along the axis of rotation at depth approximately midway between the top and bottom of 
the tank. As soon as the fluid is injected, it is deflected horizontally by the stratiflcation 
and the Coriolis force organizes it into a freely decaying anticyclone. In a second set of 
experiments, the vortex is permanently sustained by a smooth and stationary injection, 
using a peristaltic pump whose flux rate is chosen between 6 and 500 mL/min, through 
the 3 mm diameter pipe with a piece of porous material fixed at the end of the pipe. The 
technique in the second second set of experiments allows us to create vortices with higher 
Rossby numbers than in the first set. Moreover the sustained anticyclones in the second 
set of experiments evolve very slowly compared with those in the first set. In both sets 
of experiments, the injected fluid is seeded wit h fluorescein dye and 10 yitm-diameter 
particles for Particle Image Velocimetry (PIV) (JMeunier &: Lewekell2003l ) . We follow the 
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Figure 2. Image-processed side-view of a vortex boundary (white) at different times t with 
/ — 6.8 rad/s and iV = 1.6 rad/s held fixed. |_Ro| decreases in time. Also shown are the 
theoretical boundaries (gray in print; red online) of the model vortices at each time. The model 
vortex has H/L given by equation (|2.3[) . so the coincidence of the gray (or red) and white curves 
validates our law. T = 47r// is the period of the rotating turntable. Because //iV = 4.25 is the 
same for all of these vortices, while a differs, it is clear that the scaling law a = f /N is not 
correct for these vortices. 



evolution of the vortices using one horizontally and one vertically illuminated laser sheet. 
Video images of the vertical cross-section allow us to record the changes in time of the 
vortex aspect ratio (figure [IJa) and (b)), while the PIV measurements in the horizontal 
sheet allow us to find the azimuthal and radial components of the velocity of the vortex, 
which lead to the determination of Ro (figure [Ijc) and (d)). 

4.1. Freely decaying vortices 

The brief injection of fluid with density po at height z (where p{z) — po) is immediately 
followed by fast adjustments where the injected fluid becomes approximately axisymmet- 
ric. After axisymmetrization, \Vg\ and \Ro\ decay very slowly in time with the vortices 
persisting for 1000 to 1800 table rotation times (47r//), or several hundred turnaround 
times of the vortices {2tt / {f Ro)) . The vortex persists as long as the density anomaly 
{p — p) does. During the slow decay, the vortex core passes through a series of quasi- 
equilibrium states where it has approximately solid-body rotation with negligible radial 
velocity as shown on figure [TJd). The experimentally measured boundaries of one of the 
slowly-decaying anticyclones that were created by injection are shown at four different 
times in figure [21 Also shown are the boundaries of our theoretical model of the decaying 
vortex. Our model is the solid-body rotating vortex with a discontinuous velocity derived 
in section 2 with Nc = and with the ellipsoidally-shaped boundary given by equa- 
tion (12.21) . Due to the slow diffusion of salt, we assume that for all time the fluid density 
within the model vortex remains at its initial value of po, so we assume that Nc = for 
all time. Due to the lack of diffusion in the laboratory vortex, we also assume that within 
the ellipsoidal boundary of the model vortex, the volume 4/3nH{t)L{t)'^ = 4/37ra(i)L(i)'^ 
remains constant, despite the fact that a and L change in time. The theoretical bound- 
aries of our model vortices are computed with equation (|2.2[) . where a(t) is given by our 
law (|2.3p . where Ro{t) is measured experimentally, Nc is assumed to be zero, / is known, 
N is assumed to remain at its initial value, and where L{t) is computed by assuming 
that the volume of the ellipsoid is constant in time. That is, we set 



L(t)=Lo[ao/«W]'/' 



(4.1) 



where the value of Lq is the experimental value of L from the first panel with Ro = —0.08; 
ao is the initial value of a found from equation ()2.3|) with Ro — —0.08 and a{t) is 
determined from Ro{t) using equation ()2.3p . 
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Figure [2] shows that the boundaries of the laboratory and model vortices are nearly 
coincident at four different times, and at those times the vortices have four different 
Rossby numbers and four different aspect ratios. This validates the law (|2.3p for a and 
also shows that our assumptions that Nc = and that the volume of the vortex remains 
constant are both good. If the correct scaling law were H/L = f/N, rather than law (12. 3p . 
then the vortex in figure [2] would have the same a spect ratio through time, which it 
clearly does not. Using the scaling law of lGilll ( 19811 ) with a Ro dependence would lead 



to a smaller predicted aspect ratio by a factor between 3.4 and 6. 

4.2. Vortices sustained by continuous injection 
In a second set of experiments, the vortices are sustained by a continuous inie c tion o f fluid 



with density po at a fixed flux rate, as in experiments of lGriffiths fc LindenI (|l98lh . The 
characteristic values of the Ro of these vortices range between —0.45 and —0.20. With 
continuous injection the volume of these vortices slowly increases in time and the Rossby 
number of the vortices decays very slowly compared to the viscously decaying vortices 
discussed in section 4.1. The continuous injection is laminar and the rate of injection is 
sufficiently small so that the volume of the vortices increases by only approximately 1% 
per table rotation. As a consequence, we use the same theoretical, elliptical model for 
these vortices as we used for the viscously dissipating vortices analyzed in section 2. The 
experimentally measured boundaries at four different times of a vortex sustained with 
continuous injection are shown in figure [3l As before, the aspect ratio of the extracted 
shape of the laboratory vortex is compared to the a of our model vortex that uses 
law ()2.3p with Nc — 0, the observed value of Ro, for a, and the volume calculated with 
the fiux rate of the experiment and time t/T. As can be seen, the comparison is excellent 
and validates our theoretical law (j2.3p in the cyclo-geostrophic regime. 

Note that in the similar experiments of G riffiths fc Linden (^198^1) briefiy described 
in section 5 of their article, continuous intrusions of fluid into density gradients were 
performed. They observed vortices with aspect ratios of 0.47 and 0.99. In these two 
experimental runs, NH/ fL was estimated around 0.3, which is in disagreement with 
Charney's QG law. As no precise measurement of the R ossby num ber was done, no 



comparison is possible with either our law or Gill's theory (|Gilllll981l) . 



5. Application to Meddies and Jovian Vortices 

Our law (j2.3p for vortex aspect ratio a also applies to ocean meddies, which, unlike our 
laboratory vortices, are internally stratified (i.e. iVc / 0). Using the report ed values of 
the ocea n densitie s as function s of position w ithin and outsid e the med dies (^Ar mi et 



[1988; H ebert et al . 1990; Pingr ee fc Le Cann 1993; Tvchensk v fc Garto n 1998), we have 
compiled in table 1 average values of the buoyancy frequencies of the meddies Nc and 
their background environments N, along with the observed values of a, L and Ro. Note 
that according to our law (|2.3p . the effects of non-zero Nc do matter on the aspect ratio, 
especially when Nc is of order N, as it is for the meddies (and Jovian vortices - see below), 
has a large effect on the aspect ratio. For all but the oldest meddy shown in figure IH 
law (|2.3|) fits the observations. Note that setting a equals to the alternative scaling f/N 
does not fit the meddy data, even with the large uncertainties of the observed values 
of a of the meddies shown in figure HI Setting a equals to the other alternative scaling 
discussed in section 3, Ro{f /N) ([Gill 1981; He dstrom fc Armilll988l) . is an even poorer 
fit. 

Law (j2.3|) also applies to Jovian vortices. For the Great Red Spot (GRS) it is necessary 
to take into account the fact that the characteristic horizontal length L of the derivative 
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Figure 3. Image-processed side-view boundaries of a vortex sustained with continuous injection 
at different times. Tlie flow has / — f.62 rad/s, N = 2.3 rad/s, and an injection rate of 10 
mL/min. Also shown are the theoretical boundaries of the model vortices. Line colors are as in 
figure 2. The model vortex has H/L given by equation (12. 3p . so the coincidence of the boundaries 
of the laboratory and model vortex validates equation (|2.3|) for a. Unlike the vortex in figure 2, 
the vortex volume V changes in time. Ro, V, and time t are given for each image. 



of the pressure anomaly (j)—p) is nearly 3 times small er than the characterist ic radius Ry 
where the azimuthal velocity reaches it s peak v a lue ([Shettv fc Marcusll2010l ). As shown 
in our companion paper bv Hassanza deh et all ( 2012r i. when Ry ^ L, the scaling law 
should be modified so that the numerator Ro{l + Ro) in equation p.3p is replaced with 
Ro [1+Ro {L/ Ry)]. With the exceptions o f J? and iVc, the propert ies of Jovian vortices are 



well known and have small uncertainties ( Shettv fc MarcudbOlOfl . Based on observations 
of the haze layers above the GRS and Oval BA, most observers ag r ee that the elevations 



of the haze layers above the UKb and Uyal .dA, most observers ag r ee tfiat the elevations 
of their top boundaries (JBanfield et al\ Il988t iFletcher et al\ bOlOt Ide Pater et al\ 120101 ) 



are near the elevation of 140 mb pressure. There is l ess agreement on the elevation s of the 
mid-plane (z = 0) of the vortices. Some modelers ( Morales- Juberias et a/.l 120031 ) of the 
Oval BA, set z = to be 680 mb (the height of the clouds f rom which the velocit ies are 
extracted), making H — 34 km. However, some observers ( de Pater et al\ 2010 ) argue 
that z = is deeper at 2000 mb, making if = 59 km. Other modelers (jCho et aLll200lh 
choose H of the GRS to be approximately one pressure scale height (23 km). Based on 
these observations and arguments, we set the "observe d" value of H to be 45 ± 13 km. 

Jovian values of N have been measured accurately (jShettv fc Marcud 120101 ) . and as 
shown in the online supplementary material, the values of {N^ — N^) and their uncer- 
tainties for the GRS and the Oval BA can be inferred from thermal imaging, which 
provides the temperatures of the vortices and of the background atmosphere. We show 
in the online supplementary material how the values and uncertainties of Nc for the 
GRS and Oval BA that are used in table 1 are calculated from the observed tempera- 
ture measurements and also how the values of Ro for the Jovian vortices in table 1 were 
calculated. Using the values in table 1 for the properties of the Jovian vortices and of 
the background Jovian atmosphere we have included the aspect ratio information of the 
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Figure 4. Tests of our theoretical law (|2.3p (straight line) for a. The theory agrees with 
our freely decaying vortex experiments (■) with 1.6 rad/s ^ / ^ 7 rad/s and with our 
sustained vortex experiments with f = 1.6 rad/s and A^ — 2.3 rad/s (□); with Bush & 
Woods' experiments at Ro = —0.5 ("Bush & Woods" 199!f) (T); with Jupiter's Great Red Spot 
and Oval BA (Marcus 1993; Shottv & Marcus 2010; dc Pater ei al. 2010) ( A ); an d with med- 
dles Sharon (Armi et al. 1988; Hebert et al.\ 11990 : ISchlu tz Tokos & Rossby 1990)(x), Bobby 
CPing ree & Le Cann 1993) (•), Hyperion (o), Ceres (•) and Enccladc (o) (Tvchcnskv & CartoiJ 
Il99^ . All the data for the meddles are given in table 1. Data and error bars for the Jovian 
vortices are discussed in the online supplementary material. Error bars for the meddles are 
dominated by uncertainties in {N^ — N^), which are due to uncertainties in the reported den- 
sities and temperatures of meddles. Meddy Sharon is the oldest meddy shown, so dissipation 
could have modified its shape. 



GRS and Oval BA in figure |31 The figure shows that the aspect ratios of the Jovian 
anticyclones are consistent with our law (|2.3p for a (with the correction due to fact that 
Rv j^ L for the GRS). Even with the large uncertainties in H, the data show that the 
aspect ratios of the Jovian vortices are not equal to f/N or to Ro (f/N). 



6. Conclusions 

We have derived and verified a new law (j2.3|) for the aspect ratio a = H/L of vortices 
in stratified, rotating flows; the law depends on the Coriolis parameter /, the Rossby 
number Ro, the stratification of the ambient fluid N, and also on the stratification 
inside t he vortex Nc- Thi s law d erived in a more general context in the companion pa- 



per of iHassanzadeh et al\ (|2012l ) fits exactly the equilibrium solution of the Boussinesq 
equations that we used to model our laboratory anticyclones. We have shown that the 
law works well for predicting the aspect ratios of freely decaying anticyclones in the 
laboratory, of laboratory vortices that are sustained by continuous fluid injection, of 
Jovian vortices, and of Atlantic meddles. These vortices span a large range of Rossby 
numbers, occur in different fluids, have different ambient shears, Reynolds numbers and 
lifetimes, and are created and dissipated by different mechanisms. For these vortices we 
have demonstrated that a previously proposed scaling law, H/L — f/N, is not correct. 
Our law p.3p for a shows that equilibrium cyclones must have N < Nc, that anticyclones 
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iVc 


/ 


Ro 


Observed 




Data 


Vn^ - Nj" 


L 




(rad/s) 


(rad/s) 


(rad/s) 




a 


(km) 


Experiment 1 


1.6 





6.8 


-0.32 - -0.02 


1.71 - 0.68 


_ 


Experiment 2 


1.6 





3.2 


-0.16 - -0.03 


0.82 - 0.38 


- 


Experiment 3 


2.3 





1.6 


-0.45 - -0.2 


0.36 - 0.27 


- 


Bush & Wood's exp. 


0.78 - 1.4 





0.4 - 2.2 


-0.5 


0.13- 1.32 


- 


Meddy Ceres 


0.0022 


0.002 


8.4610"^ 


-0.05 


0.0073 


35 


Meddy Hyperion 


0.0024 


0.0018 


8.46 10^'' 


-0.07 


0.009 


45 


Meddy Encelade 


0.0026 


0.0015 


8.110-^ 


-0.08 


0.009 


45 


Meddy Sharon 


0.0021 


0.0011 


7.710-5 


-0.15 


0.01 


32 


Meddy Bobby 


0.0015 


0.0017 


8.3 lO-"' 


-0.17 


0.018 


27 


Jupiter's Oval BA 


0.0048 


0.018 


1.9210-'* 


-0.16 


0.014 


3200 


Jupiter's GRS 


0.0047 


0.015 


1.3710-'' 


-0.38 


0.020 


2300 
{Rv=6260 km) 



Table 1. Data for figure 3. The three sets of laboratory experiments listed in the top three 
lines have different values of N and /. Experiment 1 and 2 correspond to a freely decaying 
vortex and experiment 3 to a sustained vortex. Each set of experiments had vortices with 
different values of Ro, but in all of our laboratory experiments, we assumed Nc = and that 
Rv = L. 21/ and 2H were defined to be the horizontal and vertical diameters of the vortices 
as determined by our edge-finding algorithm (c.f. figure 2). The uncertainties in the laboratory 
measurements are dominated by the small non-axisymmetric component of the vortices. The 
fractional uncertainties in iV are ±10%. For the meddles, the values of A'' and Nc are such that 
1.33 < N /Nc < 2.16. The parameter values and their uncertainties for the Jovian vortices are 
discussed in section 5 and in the online supplementary material. Note that the Jovian values of 
N, f and L are known with small uncertainties but that H has large uncertainties. 



with \Ro\ < 1 have N > Nc and that anticyclones with \Ro\ > 1 have N < Nc- Mixing 
within a vortex tends to de-stratify the fluid inside of it and therefore naturally decreases 
Nc from its initial value. Therefore, if there is mixing and if there is no active process that 
continuously rc-stratifies the fluid within a vortex, then at late times one would expect 
that Nc < N, and then according to our law only anticyclones with Ro < 1 can be in 
equilibrium. This may explain why there are more anticyclonic than cyclonic eddies in 
the ocean, and also why there appear to be more long-lived anticyclones than cyclones 
on Jupiter, Saturn and Neptune. 
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7. Properties of the Jovian atmosphere and vortices 

To compute the values of the GRS and the O val BA shown in table 1, we used g 



25 m s"2 and iJ = 45 ± 13 km. For the GRS (|de Pater et all l2010l : IShettv fc MarcusI 
20101) we used N = 0.0158 ± 0.0005 rad/s, Ve = 119 ± 14 m/s, L = 2300 ± 70 km, 
Ry = 6260 ± 50 km, and / = 1.374 x lO"" rad/s. For the Oval BA we used N = 
0.0182 ± 0.0006 rad/s, Vf, = 101 ± 9 m/s, L = 3200 ± 45 km, R„ = 3200 ± 45 km, and 



/ = 1.915 X lO-'' rad/s (Ide Pater et aLll2010l : [Shettv fc Marcusll2010l) . 

Thermal imaging of the Jovian atmosphere gives the values of the temperature of the 
clouds at the elevations of the Great Red Spot (GRS) and Oval BA. In particular, the 
imaging gives the temperature Tc' (where c' indicates the intersection of the vertical 
central rotation axis of a vortex and the the top of the vortex at z = H). Thermal 
imaging also gives the value of the background temperature T\z=h of the atmosphere at 
the elevation of the top of the vortex and at the latitude of the principal east-west axis 
of the vortex. For a vortex that is in cyclo-geostrophic equilibrium, at all elevations z, 
the Coriolis and centrifugal accelerations due to the azimuthal velocity of the vortex ar e 
balanced by the horizontal pressure gradient within the vortex ( Hassanzadeh et aLll2012r ). 
The top of the vortex, z = H, is defined to be the location where the azimuthal velocity 
of the vortex is zero. Therefore at z = _ff , the horizontal pressure gradient within the 
vortex is zero. Therefore, in an ideal gas at the top of a vortex, [Tc' — {T\z=h)]/{T\z=h) = 
-[pc'-ipU=H)]/ipU=H) ^ -H[idp/dz)U'-idp/dz)]/ipU^H) ^ iH/g)iN^-N^), where 
the first approximate equality in this expression comes from using the first term in a 
Taylor series expansion to compute the density difference, and the second approximate 
equality comes from the definitions of the buoyancy frequencies. Thus, we can set the 
value of (N^ — N^) for the Jovian vortices equal to the observed value of ^* ^Z}, '^"!^'' . 

Following [Morales- Juberias et al\ (|2003() , we set the upper elevation of the GRS to be 



at 140 mbar and use temp erature measuremen ts at that elevation. Satellite observations 
of the GRS (see figure 8 of lFletcher eraLll2010l) give [T^' - iT\z=H)]/{T\z=H) = -0.041 ± 
0.015. Our estimate of the uncertainty of the Jovian temperatures of the GRS are based 
on the standard deviation of the data in the top-right panel of figure 8 in iFletcher et all 
(|2010[ ). The values of ^/W^^~^W^ and N^ in table 1 for the GRS are based on these 



temperatures and uncertainties. 

There have been no direct measurements of the temperature difference across the top 
of the Oval BA, but there were thermal measurements of the three White Oval vortices 
that existed at the same latitude as the Oval BA before it formed (and from which it 
formed) . We use thes e White Ovals as a proxy for the Oval BA because it has been argued 
( de Pater et aZ.ll2010l ) that the White Ovals were dynamically similar to the Oval BA and 
also of the same size and shape. For the Ovals, [T^' - {f\,,= H) ]l{f\r.= H) = -0.04 ±0.015. 



For the Oval BA, we used the temperatures from figure 1 of lConrath et al\ (|l98lf ). The 
uncertainty of those temperature measurements was not published, so our estimate of 
the uncertainty is based on the standard deviation of the values of \Tc' — {T\z=h)\ of 
the three White Ovals at 140 mbar. The values of \/N'^ - N^^ and Nc in table 1 for the 
Oval BA are based on these temperatures and uncertainties. 

To apply our equation (J2.3I) correctly to the Jovian vortices, which are not axisymmet- 
ric because they are embedded in strongly shearing east-west winds, it is necessary to 
consider the derivatio n of the equation for a bas ed on cyclo-geostrophic balance shown in 
the companion paper (JHassanzadeh et alll2012| ). The law for a can be obtained from the 
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vertical hydrostatic equilibrium and cyclo-geostrophic balance applied to the east-west 
component of the force along the principal east-west axis of the vortex. The shearing east- 
west wind does not enter into this balance and therefore the Rossby number Ro — VgL/ f 
that comes from this derivation and that appears in equation (I2.3P for a (and that is re- 
ported in table 1), is based on the values of the north-south velocities along the east-west 
principal axis. Similarly the value of L that should be used in Ro (and that is reported 
in table 1) is that of the characteristic length scale of the pressure gradient along the 
east- west principal axis. Consistent with our approach of deriving the relationship for a 
by balancing the east-west component of the forces along the east-west principal axis, 
the measurements of [Tc' — {T\z=h)] that should be used in approximating {N^ — N'^) of 
the Jovian vortices were made along the vortex's east-west principal axis. 



